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Abstract 

In this paper, the Bosonic projected variational method is proposed to 
study the strongly correlated 87 Rb atoms in optical lattice. A global phase 
diagram is obtained by this method. There exist two characteristic lattice 
depths V for 87 Rb atoms in optical lattice : one is V = 9.5E r to label the 
maximum height of the 'zero-momentum' peak of condensation, the other 
is the quantum critical point for the superfluid-insulator (SI) transition at 
V = 12.3£y. As a result of strongly correlated effect for lattice Bosons, 
the suppressed superfluid state is predicted near the SI transition with the 
suppressed superfluid density and the very slowly velocity of the sound-like 
excitons. 
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During the last years a spectacular development in the storage and manipulation of 
cold atoms in optical lattices has taken place. An superfluid-insulator (SI) transition was 
observed for 87 Rb atoms trapped in a three-dimensional optical lattice potential by changing 
the potential depth [1]. The SI transition is the appearance of a excitation gap. In the 
superfluid regime, there is no excitation gap and instead one observes a series of 'Bragg- 
peaks' around the characteristic 'zero-momentum' peak of a condensate in the absence of 
an optical lattice. The existence of this gap has been verified experimentally by applying a 
phase gradient in the Mott-insulator. This experimental progress has revived the interest in 
the Bose- Hubbard (BH) model as a generic Hamiltonian for strongly correlated bosons [2], 
by which the quantum phase transition can be described [3-18]. 

In this paper we intend to examine the physics in 87 Rb atoms by studying a Bose-Hubbard 
model using a variational method for bosonic systems [19]. In the new approach, the on-site 
repulsion is treated exactly, while the kinetic energy is studied variationally, so that it is 
suitable to examine some issues in strongly correlated Bosonic systems. The variational 
method applied to the Bose-Hubbard model in three-dimension demonstrates a quantum 
phase transition from a superfluid with suppressed superfluid density for smaller intra-site 
Coulomb repulsion U to a Mott insulator for larger U at unit filling ^ = 1 . There exists 
a quantum critical point (QCP) for the homogenous phase. We calculate the superfluid 
density near the QCP and show the existence of the suppressed superfluid state which is a 
kind of "gossamer" phenomenon. The idea of "gossamer" (superconducting) state is that the 
"insulator" might actually be a thin, ghostly superconductor which is proposed by Laughlin 
[20,21]. In a suppressed superfluid state, the superfluid density is very thin and the velocity 
of the sound-like excitons is very slow, in contrast to the conventional superfluid state. 

The system is based on confining cold 87 Rb atoms in the periodic potential of an op- 
tical lattice [1]. In the simplest case, three orthogonal, independent standing laser fields 
with wave vector k produce a separable three dimensional lattice potential V(x,y,z) = 
V (sin 2 kx + sin 2 ky + sin 2 kz^j with a tunable amplitude V ^> E r = h 2 k 2 /2m. Starting 
from the standard pseudopotential description the interatomic potential is replaced by an 
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effective contact interaction of the form U(x) = 47r ^° s • 5(x) containing the exact s-wave 
scattering length a s . With a\ as the creation operator of a boson at site i and hi as the 
density operator, the Hamiltonian reads 

H = - J h % + J 12 - X ) (!) 

<ij> i 

The bandwidth parameter J is essentially the gain in kinetic energy due to nearest neighbor 

A ( V \ 3 / 4 

tunneling J = -^E r ^ J-J exp —2 ^-^-J . The relevant interaction parameter [/ is thus 
given by an integral over the on-site wave function via U = ^J~^ka s E r (Jf) ^ ■ 

We consider \^>) = P p \ $ ) as a variational wavefunction to examine the ground state 
of the above effective BH model for 87 Rb atoms in optical lattices. | $o) is the wavefunction 
of the Bose-condensed ground state for non-interacting Bosons 

|$o >= exp(v / iVe^°aj {=o )|0 > (2) 



where iV is the number of Bosons. The order parameter ($o|^ =0 |$o) = e m WN shows an 
off diagonal long range order. P p is the Bosonic partial projection (BPP) operator defined 
as 



N N 



m=2 i j=0,j^m V' 1 J) 

If the on-site repulsive interaction between bosons is large enough, there exist little possibly 
of high occupation states, we can introduce the Bosonic partial projection (BPP) operator 
to describe the strongly correlated effects. If e rn = 0, the Bosonic partial projection operator 
P p turns into the Bosonic completely projection operator P c , the ground state | ^) is reduced 
into a Mott insulator state with gapped excitons. If e m — 1, there is no interaction between 
Bosons, the Bosonic partial projection operator P p turns into constant number P p — 1, and 
the ground state turns into a Bose-condensation state with massless excitons E ~ aq 2 . 
The variational energy for the ground state becomes 

E g = (* | H | | *) (4) 

= (H t ) + UD 2 + 3UD 3 + ... + iV(iV ~ l) UD N . 



The occupation numbers are D 2 , D 3 , ...D N . The variational kinetic energy (H t ) is 

(H t ) = (* | -JE(«1% + ^-c) I I = 9jE (5) 

(ij) 

where E the kinetic energy for non-interaction Bose systems 

E = ($ | -J£(«!% + ^.c.) | $ )/($o I $o)- (6) 

gj is renormalization factor for kinetic energy to be determined. 

In the thermodynamics limit iV — > oo, we choose the maximum term of the variational 
wavefunction | ^) and obtain the condition 

2 = d m (x + d 2 + 2d 3 + ... + (N — l)^)"- 1 
£m (l-x-2rf 2 -3rf 3 --^ i v) m ' 

m = 2, 3, ...N. 

where x = (N — D 2 — D 3 — ...D N )/L describes the empty concentration. L is the number 
of the site. From the condition we have the renormalization factor for the kinetic energy 

N j^m 

9J=J2 C m-^Z^, (8) 
m=l D 

2^m-l + ^2^3^-3^-2 + ••• 
~^~^i^m—i^i—l^rn—i+l ■■■ ^-m—\^-m- 

with A m = 1 - x - 2d 2 - 3d 3 ... - iVdjv and5 = i + d 2 + 2d 3 + ... + (N — l)d N . 

In the strong interaction limit, one has e 2 m ~ d ^fx-2d 2 ) m ' '^' ie renorma li za tion factor for 
the kinetic energy is obtained as 

gj ~ (d 2 + x)(l - x - 2d 2 ) (9) 
+...2yjd rn -2d m -i[(d 2 + x)(l -x- 2d 2 )] 1/2 + .... 

where d m = -^f- are determined by the equations 

dE g (d m ) 



dd m 



0. 
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In this paper we solve the equations by the approximation up to d 2 without considering 
d m , m > 2 [22]. The renormalization factor is obtained as gj (1 — x — 2d 2 )(2d 2 + 



x + 2y d 2 [d 2 + x)). d 2 is determined by the minimum energy condition = as d 2 ~ 
\[l - In the following paper we consider a BH model with unit filling N/L = 1 or 

x = 0. 

Note that in the BPP approach and at unit filling case N/L = 1 or x = 0, d 2 is a measure 
of the mobile carrier density. At d 2 = 0, we have (H t ) = 0. This state describes a Mott 
insulator. On the other hand, the case with d 2 > describes superfluid state. We expect a 
transition from the Mott insulator at larger U to the superfluid at smaller U as U decreases 
passing through a critical point U c . The transition point U c is given by 

*.= <^# = *7 (10) 

or (^-) = i ln 2 (-^|p) ~ 12.3. For U > U c = 24J, there is no solution for physical values of 
d 2 , indicating that d 2 = 0. 

If U < U c ( (^) < 12.3), d 2 > (e 2 > 0), the ground state is a superfluid state. 

Within the BPP approximation, the suppressed superfluid density ($|ng =0 |$) is 

m = i[i - (ID 

In this state, the mobile carrier density d 2 is fixed and the excitons are the phase fluctuations 
- phasons. To consider the dynamics of phasons, the quantum states are defined as 

L 

|* >= P p expQ>^a+)|0 > ■ ( 12 ) 

3 

The effective Hamiltonian becomes 

Heff = (* | -JE(«1% + fc-c) + ^ E"i("i " !) I I *> ( 13 ) 

where p p = |[1 — (^-) 2 ] is the phase stiffness. As a result the effective model of the SF state 
turns into a quantum four dimensional XY model with long range order. The transition 
temperature T c is scaled as T c ~ p 2 / 3 [2]. 
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In the SF regime, there is no excitation gap and instead the homogeneous system exhibits 
a sound like mode - phason with frequency u(q) = cq. The associated sound velocity is c = 
^JjJU[l — (jTf-) 2 ]. It is the correlated effect drives the BEC-like q 2 spectrum to superfluid- 
like q type. 

The SF phase of 87 Rb atoms in an optical lattice can thus quite generally be characterized 
by the fact that at reciprocal lattice vectors q — 0, the momentum distribution n has a 
peak with the height as [1,24] 

n = L-($|n A ; =0 |$)K0)| 2 (14) 
~ L\[l - V2ka s eM-^V/E r )](^-f 2 . 

The fact that the peaks in the momentum distribution at q = initially grow with increasing 
depth of the lattice potential is a result of the strong decrease in spatial extent of the Wannier 
function w(x), which entails a corresponding increase in its Fourier transform w(q) oc (J-) 3 ^ 4 - 
Beyond a critical lattice depth around V = 9.5E r , this trend is reversed, however, and the 
superfluid density eventually disappear completely at U — U c (V — 12.3E r ). 

If U > U c or ( ¥-) > 12.3, d 2 = and e 2 = ( d m = and e m = 0, m > 2), the ground 
state is a Mott insulator, of which the wave function is 

|$>=P c |$ >=L- 1 / 2 n(|l Z >). (15) 

I 

The energy gap describes the energy difference between the ground state with D 2 = 
and the excited state with D 2 = 1 as AE ~ — 24Jcos(ag/2) + U. The excited energy AE 
is from (U — 24 J) to (U + 24J) with the center at U. A is the Mott gap for Boson exciton 
defined as A = (U — 24 J). Deep in the MI phase, this gap has size U, which is just the 
increase in energy if an atom tunnels to an already occupied adjacent site. The existence 
of this gap has been verified experimentally by applying a phase gradient in the MI and 
measuring the resulting excitations produced in the SF at smaller V/E r [1]. Near the Mott 
transition, the gap closes A — (U — 24J) = 0. The signal of the Mott gap A disappeared 
around V = 12.3E r (U = 24 J), which was taken also as another definition of the critical 
point of the SI transition. 
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Above results predict two characteristic lattice depths: V = 9.5E r and V = 12.3E r . 
Near the lattice depth around V = 9.5E r , the 'zero-momentum' peak has maximum height 
at which the condensation is most robust. Beyond this point, the phase coherence becomes 
weak. The other character lattice depth is V = 12.3E r which is just the QCP for the SI 
transition. At this point the 'zero-momentum' peak disappears together with the opening 
of a gap for particle-hole excitons. From the experiments by M. Greiner, at almost the same 
point V ~ 13E r the 'zero- momentum' peak of a condensate has maximum height and the 
Mott gap A opens. Thus there exist only one characteristic lattice depth. The difference 
between our results and the experiments is due to the inhomogeneous [1,4,25]. In this paper 
our theory is based on the homogenous phase without considering the inhomogeneous phase. 
The inhomogeneous phase in the BH model from BPP method will be explored elsewhere. 

In strongly correlated limit near the SI critical point < U u Uc C 1, a new class of 
superfluid state appears - the suppressed superffuid state. The wavefunction for so called 
the suppressed superfluid state is read as 

|$ > = P p \% >= L- 1 / 2 JI(|0i >+|l*> + E >), (16) 

i n l= i n i- 

d m {x + d 2 ) m ~ 1 

F ~ <ST 

ni ~ (1 - x - 2d 2 )™ ^ 
It is obviously that the superfluid density for the suppressed superfluid is suppressed 
seriously by correlations 

n(q = 0) = ($| Vo|$) = ^(1 " (^) 2 ) - 0, (17) 

a quantitative measure of the suppressed superfluid. Another feature for the suppressed 
superfluid is the existence of a very slowly velocity of the phasons. The associated sound 
velocity turns into zero c = ^J^JU[1 — (jy-) 2 } — > 0. The third character for the suppressed 
superfluid is the pinned chemical potential at the center of the gap. 

Let us now consider the evolution of chemical potential. The chemical potential /i is /i — 
+ Tfcf ■ The term originates from the x dependences of gj in the variational procedure [21], 
which will be important in calculation of the chemical potential of the state 
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dgj 



1 - 2(2d 2 +x)- 2^d 2 (d 2 + x) 



(18) 



dx 



+ (1 - 2d 2 - x)^{d 2 + x)/d 2 . 



In the limit x 



8^2, one has 



U 



(19) 



Our results show that the suppressed superfluid state is a similar phenomenon to the 
gossamer superconductivity. From the results of the partial projection to the Gutzwiller 
variational method, the ground state for strongly correlated electrons may be superconduct- 
ing at half filling due to some kinds of attraction mechanism [21]. In this paper it is shown 
that because of the partial projection, there exists the suppressed superfluid state at unit 
filling j; — 1 in the region of < u ~ fJ Uc <C 1. The Bose condensation can be suppressed 
seriously by strongly correlated effect. 

In summary, we have used the Bosonic projected variational method to study the Bose- 
Hubbard model which is the effective model for strongly correlated 87 Rb atoms in the optical 
lattice. And we obtain a global phase diagram of the homogenous phase and give two results 
: two characteristic lattice depths in phase diagram and the existence of the "gossamer" -like 
state near the SI transition. 



The author thanks W.M. Liu for helpful conversations. Kou Su-Peng. acknowledges 
support from NSFC Grant No. 10574014 and Beijing Normal university. 



ACKNOWLEDGMENTS 



8 



References 



[1] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and I. Bloch, Nature 415, 39 (2002). 

[2] M.P.A. Fisher, P.B. Weichman, G. Grinstein, and D.S. Fisher, Phys. Rev. B 40, 546 (1989). 

[3] S. Sachdev, Quantum Phase Transitions, (Cambridge Univ. Press, Cambridge, 2001). 

[4] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P. Zoller, Phys. Rev. Lett. 81, 3108 
(1998); D. Jaksch, V. Venturi, J. I. Cirac, C. J. Williams, and P. Zoller Phys. Rev. Lett. 89, 
040402 (2002). 

[5] C. Bruder, R. Fazio, and G. Schon, Phys. Rev. B 47, 342 (1993). A.P. Kampf and G. T. 
Zimanyi, Phys. Rev. B 47 279 (1993). K. 

[6] Sheshadri, H.R. Krishnamurty, R. Pandit, and T.V. Ramakrishnan, Europhys. Lett. 22, 257 
(1992). 

[7] D.S. Rokhsar, B.G. Kotliar, Phys. Rev. B 44, 10328 (1991). 

[8] W. Krauth, N. Trivedi, and D. Ceperley, Phys. Rev. Lett. 67, 2307 (1991). W. Krauth and 
N. Trivedi, Europys. Lett. 14, 627 (1991); W. Krauth, M. Caffarel, and J.-P. Bouchard, Phys. 
Rev. B 45, 3137 (1992). 

[9] R.T. Scalettar, G.G. Batrouni and G.T. Zimanyi, Phys. Rev. Lett. 66 3144 (1991) 

[10] P. Niyaz, R.T. Scalettar, C.Y. Fong and G.G. Batrouni, Phys. Rev. B 44, 7143 (1991) 

[11] G.G. Batrouni and R.T. Scalettar, Phys. Rev. B 46, 9051 (1992). 

[12] J.K. Freericks and H. Monien, Europhys. Lett. 26, 545 (1993). 

[13] W. Zwerger, Journal of Optics B 5, S9 (2003). 

[14] K. Goral, L. Santos and M. Lewenstein, Phys. Rev. Lett. 88, 170406 (2002). 

9 



[15] D. van Oosten, P. van der Straten, H.T.C. Stoof, Phys. Rev. A 63, 053601 (2001). D.B.M. 
Dickerscheid, D. van Oosten, P.J.H. Denteneer and H.T.C. Stoof, Phys. Rev. A 68, 043623 
(2003). 

[16] L. Amico and V. Penna, E-print: cond-mat/9801074. 
[17] N. Elstner and H. Monien, E-print: cond-mat/9905367. 

[18] Jun-Jun Liang, J.-Q. Liang, and W.-M. Liu, Phys. Rev. A 68, 043605 (2003). 
[19] M.C. Gutzwiller, Phys. Rev. A 137, 1726 (1965). 
[20] R.B. Laughlin, E-print: cond-mat/0209269. 
[21] F.C. Zhang, Phys. Rev. Lett 90, 23 (2003). 

[22] In strong interaction limit U » J, c?3 is much smaller than d,2, ds oc f 0^(1 — 2g?2) << cfe- 
[23] V.A. Kashurnikov, N.V. Prokof'ev, B.V. Svistunov, Phys. Rev. A 66 031601 (R) 2002. 
[24] R. Roth and K. Burnett 2002, Preprint cond-mat/0209066. 
[25] G.G. Batrouni et al., Phys. Rev. Lett, 89, 117203 (2002). 

Figure Captions 

Fig.l 

This figure shows the double occupation rate d 2 and the triple occupation rate d 3 of 

I E - 

Fig.2 

Phase diagram for the Bose-Hubbard model under unit filling. In this figure the scales at 
left axis and right axis are T C (U — > 0) and U c , respectively. -j^- is the parameter for Bosons 
in optical lattice. 



10 



This figure "becfigl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/05 1 2008 v 1 



This figure "becfig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/05 1 2008 v 1 



